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Abstract. An approach to compatibility analysis of systems of discrete 
relations is proposed. Unlike the Grobner basis technique, the proposed 
scheme is not based on the polynomial ring structure. It uses more prim- 
itive set-theoretic and topological concepts and constructions. We illus- 
trate the approach by application to some two-state cellular automata. 
In the two-state case the Grobner basis method is also applicable, and 
we compare both approaches. 

1 Introduction 

A typical example of a system of discrete relations is a cellular automaton. 
Cellular automata are used successfully in a large number of applications. 1 Fur- 
thermore, the concept of cellular automaton can be generalized, and we consider 
the following extension of the standard notion of a cellular automaton: 

1. Instead of regular uniform lattice representing the space and time in a 
cellular automaton, we consider more general abstract simplicial complex 
K = (X,A) (see, e.g., Here X = {xq,x\, . . .} is a finite (or countably 
infinite) set of points; A is a collection of subsets of X such that (a) for all 
Xi £ X, {xi} e A; (b) if r C S E A, then r E A. 

The sets {x{\ are called vertices. We say 5 6 A is a k— simplex of dimension 
k if \5\ = k + 1, i.e., dim S = \5\ — 1. The dimension of complex K is defined 

as the maximum dimension of its constituent simplices dim K = max dim 5. 

SeA 

If r C S, t is called a face of S. Since any face of a simplex is also a sim- 
plex, the topological structure of the complex K, i.e., the set A is uniquely 
determined by the set of maximal simplices under inclusion. 
One of the advantages of simplicial complexes over regular lattices is their 
applicability to models with dynamically emerging and evolving rather than 
pre-existing space-time structure. 

1 Comparing expressiveness of cellular automata and differential equations, T. Toffoli 
writes [JJ : "Today, it is clear that we can do all that differential equations can do, and 
more, because it is differential equations that are the poor man's cellular automata 
— not the other way around!" 



2. The dynamics of a cellular automaton is determined by a local rule 

x ik = / \ x io> ■ ■ ■ i x ik-i) ■ (1) 

In this formula Xi , . . . , Xi k 6 X are interpreted as discrete variables taking 
values in a finite set of states S canonically represented as 

S = {0,. ..,<?-!}. 

The set of points {xi , . . . ,Xi k _ 1 } is called the neighborhood. The point Xi k 
is considered as the "next time step" match of some point, say Xi k _ 1 , from 
the neighborhood. 

A natural generalization is to replace function by a relation on the set 
{xi , . . . ,Xi k } . In this context, local rule JIJ is a special case of relation. 
Relations like Q ar e called functional relations. They are too restrictive in 
many applications. In particular, they violate in most cases the symmetry 
among points Xi , . . . , Xi k . Furthermore, we will see below that the functional 
relations, as a rule, have non-functional consequences. 

We can formulate some natural problems concerning the above structures: 

1. Construction of consequences. Given a relation R s on a set of points 8, 
construct non-trivial relations R T on subsets r C 8, such that R s =4> R T . 

2. Extension of relation. Given a relation R T on a subset r C 8, extend it to 
relation R s on the superset 8. 

3. Decomposition of relation. Given a relation R on a set 8, decompose R into 
combination of relations on subsets of <5. 

4. Compatibility problem. Given a collection of relations {R Sl , . . . , R s " } defined 
on sets {Si, . . . ,8 n }, construct relation R u ?=i Si on the union U"=i sucn 
that R u i=i Si is compatible with the initial relations. 

5. Imposing topological structure. Given a relation R x on a set X , endow X 
with a structure of simplicial complex consistent with the decomposition of 
the relation. 

If the number of states is a power of a prime, i.e., q = p n , we can always 2 
represent any relation over k points {xi, . . . ,Xk} by the set of zeros of some 
polynomial from the ring ¥ q [xi,... ,Xk] and study the compatibility problem 
by the standard Grobner basis methods. It would be instructive to look at the 
compatibility problem from the set-theoretic point of view cleared of the ring 
structure influence. 

An example from fundamental physics is the holographic principle proposed 
by G. 't Hooft and developed by many authors (see |4I5|). According to 't Hooft 
the combination of quantum mechanics and gravity implies that the world at the 
Planck scale can be described by a three-dimensional discrete lattice theory with 
a spacing of the Planck length order. Moreover, a full description of events on 

2 Due to the functional completeness of polynomials over F 9 (see 0) any function 
mapping k elements of ¥ q into ¥ q can be realized by a polynomial. 
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the three-dimensional lattice can be derived from a set of Boolean data (one bit 
per Planck area) on a two-dimensional lattice at the spatial (evolving with time) 
boundaries of the world. The transfer of data from two to three dimensions is 
performed in accordance with some local relations (constraints or laws) defined 
on plaquettes of the lattice. Since the data on points of the three-dimensional 
lattice are overdetermined, the control of compatibility of relations is necessary. 
Large number of constraints compared to the freedom one has in constructing 
models is one of the reasons why no completely consistent mathematical models 
describing physics at the Planck scale have been found so far. 

2 Basic Definitions and Constructions 

The definition of abstract /c-simplex as a set of k + 1 points is motivated by the 
fact that k + 1 points generically embedded in Euclidean space of sufficiently high 
dimension determine fc-dimensional convex polyhedron. The abstract combina- 
torial topology only cares about how the simplices are connected, and not how 
they can be placed within whatever spaces. 3 We need to consider also fc-point 
sets which we call k-sets. Notice that fc-sets may or may not be (k — l)-simplices. 

A relation is defined as a subset of a Cartesian product S x • • • x S of the set 
of states. Dealing with the system of relations determined over different sets of 
points we should indicate the correspondence between points and dimensions of 
the hypercube S x • • • x S. The notation S^ x *^ specifies the set S as a set of values 
for the point Xi. For the fc-set S = {x\, . . . , Xk} we denote S 5 = S^ 1 ^ x • • -xS^ Xk ^. 

A relation R s over a fc-set S = {x\, . . . , x^} is any subset of the hypercube 
S 5 , i.e., R s C S 5 . We call the set S domain of the relation R 5 . The relations 5 
and S s are called empty and trivial, respectively. 

Given a set of points S, its subset r C 5 and relation R T over the subset r, 
we define extension of R7 as the relation 

R 5 = R T x S sv . 

The procedure of extension allows one to extend relations R Sl , . . . , R Sm defined 
on different domains to the common domain, i.e., the union 5\ U • • • U 6 m . 

Now we can construct the compatibility condition of the system of rela- 
tions R Sl , . . . , R Sm . Naturally this is intersection of extensions of the relations 
to the common domain 

rs = n { R5i x sSKSz ) ' where 6 = u 6i - 
i=i ' i=i 

We call the compatibility condition R s the base relation of the system of 
relations R Sl , . . . , R 5m . If the base relation is empty, the relations R Sl , . . . , R 5m 
are incompatible. Note that in the case q = p n the compatibility condition can 

3 There are mathematical structures of non-geometric origin, like hypergraphs or block 
designs, closely related conceptually to the abstract simplicial complexes. 
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be represented by a single polynomial, in contrast to the Grobner basis approach 
(of course, the main aim of the Grobner basis computation — construction of 
basis of polynomial ideal is out of the question) . 

A relation Q s is a consequence of relation R s , if R s C Q s C S s , i.e., Q s is 
any superset of R s . Any relation can be represented in many ways by intersec- 
tions of different sets of its consequences: 

R s = Q Tl n ■ ■• nQ r '\ 

We call such representations decompositions. 

In the polynomial case q = p n , any possible Grobner basis of polynomials 
representing the relations R Sl R Sm corresponds to some decomposition of 
the base relation R s of the system R Sl , . . . , R Sm . However, the decomposition 
implied by a Grobner basis may look accidental from our point of view and if 
q 7^ p 11 such decomposition is impossible at all. 

The total number of all consequences (including R s itself and the trivial 
relation S s ) is, obviously, 

2 (<r*-|fl { l), 

In our context it is natural to distinguish the consequences which are reduced 
to relations over smaller sets of points. 

A nontrivial relation Q r is called proper consequence of relation R s if r is 
a proper subset of 5, i.e., r C 5, and relation Q T x S S \ T is consequence of R s . 

There are relations without proper consequences and these relations are most 
fundamental for a given number of points k. We call such relations prime. 

If relation R s has proper consequences R Sl , . . . , R Sm we can construct its 
canonical decomposition 

R s = PR s f)(j\(R s < xS^)j , (2) 

where the factor PR S , which we call the principal factor, is defined as 

PR S = R s |J (s s \ fl (R s > x S s \ 5 >)j . 

The principal factor is the relation of maximum "freedom", i.e., closest to the 
trivial relation but sufficient to restore R s in combination with the proper con- 
sequences. 

If the principal factor in canonical decomposition @ is trivial, then R s can 
be fully reduced to relations over smaller sets of points. We call a relation R s 
reducible, if it can be represented in the form 

m 

R s = f| (i? 5 * x S s ^) , 
i=i 
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where all R Si are proper consequences of R s . For brevity we will omit the trivial 
multipliers in intersections and write in the subsequent sections expressions like 
f}? =1 R 5 ' instead of f}?=i (R 5i x S s ^ ) . 

We see how to impose the structure of simplicial complex on an amorphous 
set of points X — {xo,xi, . . .} via a relation R x . The maximal simplices of A 
must correspond to the irreducible components of the relation R x . Now we can 
evolve — starting only with a set of points and a relation on it (in fact, we 
simply identify dimensions of the relation with the points) — the standard tools 
of the algebraic topology like homology cohomology, etc. 

We wrote a program in C implementing the above constructions and manip- 
ulations with them. Below we illustrate application of the program to analysis 
of Conway's Game of Life [H] and some of the Wolfram's elementary cellular 
automata jx]. 

A few words are needed about computer implementation of relations. To 
specify a fc-ary relation R k we should mark its points within the fc-dimensional 
hypercube S , i.e., define a characteristic function \ '■ S — > {0, 1} , with x( s ) = 
1 or according as s 6 R k or s ^ R k . Here s = (s , s\, . . . , Sk-i) is a point 
of the hypercube. The simplest way to implement the characteristic function is 
to enumerate all the q k hypercube points in some standard, e.g., lexicographic 
order: 
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Then the relation can be represented by a string of q k bits. We call this string 
bit table of relation. Symbolically Bit Table [ i or d] '■= (s G R k ) ■ Note that s is a 
("little-endian") representation of the number i orc i in the base q. Most manipula- 
tions with relations are reduced to very efficient bitwise computer commands. Of 
course, symmetric or sparse (or, vice versa, dense) relations can be represented 
in a more economical way, but these are technical details of implementation. 



3 Conway's Game of Life 

The local rule of the cellular automaton Life is defined over the 10-set 8 
{x , ■ ■ - ,xg}: 

X 5 ) (X4) 




@) @ (g) 

Here the point Xg is the next time step of the point Xg. The state set S is {0, 1}. 
The local rule can be represented as a relation R s Li f e on the 10-dimensional 
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hypercube S s . By definition, the hypercube element belongs to the relation of 
the automaton Life, i.e., (xo, • • • , xg) G R S nf e , in the following cases: 

1- fel=o^ = 3) A (39 = 1), 



2. 



(ELo x * = 2 ) A ( x 8 = x 9 ), 



3. xg = 0, if none of the above conditions holds. 



The number of elements of R Li f e is 



E><5 

n Life 



= 512. The relation R s Li f el as is 



the case for any cellular automaton, is functional: the state of xg is uniquely 
determined by the states of other points. The state set S = {0, 1} can be ad- 
ditionally endowed with the structure of the field F2. We accompany the below 
analysis of the structure of R^ e by description in terms of polynomials from 
F 2 [xo, ■ • • , xg] . This is done only for illustrative purposes and for comparison 
with the Grobner basis method. In fact, we transform the relations to polyno- 
mials only for output. This is done by computationally very cheap Lagrange 
interpolation generalized to the multivariate case. In the case q = 2, the poly- 
nomial which set of zeros corresponds to a relation is constructed uniquely. If 
q = p n > 2, there is a freedom in the choice of nonzero values of constructed 
polynomial, and the same relation can be represented by many polynomials. 
The polynomial representing R s Li f e takes the form 

Phife = Xg + X 8 {07 + <7 6 + O3 + a 2 } + a 7 + O3, (3) 

where er^ = <7k (x , . . . , X7) is the fcth elementary symmetric "polynomial defined 
for n variables x , . . . , x„_i by the formula: 

0~k (xg, . . . , X n _l) ^ ^ X^Xjj • • • Xj fc _ x . 

0<io<ii<-<U-i<n 

The relation R Li f e is reducible. It decomposes into two equivalence classes 
(with respect to the permutations of the points xo, . . . , x?) of relations defined 
over 9 points: 

1. Eight relations i?f {a:i} , < i < 7. 

Their polynomials PI (xo, ■ ■ ■ , Xi, . . . , X7, x 8 , x g ) take the form 

Pi = x 8 x 9 {erg + <7g + a\ + cr\}+x g {erg + a\ + l}+x 8 {a\ + a\ + cr l 3 + a l 2 } , 



°~k =<Tk(XQ,...,Xi,...,X 7 ). (4) 

2. One relation i?^^ 8 ^ w ith polynomial P 2 (xq, . . . ,xr,xg): 

P 2 = Xg {(T 7 + <7 6 + (T 3 + (T 2 + 1} + 0"7 + CT 3 , (7fe EE CTfc (x , . . . , X 7 ) . (5) 
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The relation R s Li f e has the following decomposition 




where (io, . . . , i&) are any 7 different indices from the set (0, . . . , 7). 

We see that the rule of Life is defined on 8-dimensional space-time simplices. 
Of course, this interpretation is based on the concepts of the abstract combina- 
torial topology and differs from the native interpretation of the game of Life as 
a (2+l)-dimensional lattice structure. 

The relations R S ^ X '^ and R^ Xs ^ are irreducible but not prime, i.e., they 
have proper consequences. 

The relation R x has two classes of 7-dimensional consequences: 

1. Seven relations R^^ x% ' x ^ with polynomials 

( x 0> ■ ■ ■ i x ii ■ ■ ■ ,Xj, . . . , X7, X%, Xg) = 

x$x 9 |fJg J + (Trj + a l l + a l 3 j + a l J + cr[ J + 1 j 

+ X9 {a$ + a i >+a*+a%+v? + l}, (7) 
a% = a k (x , ...,£i,...,£j,...,xr). 

2. One relation R^^"* 8 ^ with polynomial 

PI 2 (X Q , . . . , Xi, . . . , X 7 , Xg) = Xg + a\ + o\ + <j| + u\ + (l\ + l} . (8) 

The 8-dimensional relation R^ Xs ^ has one class of 7-dimensional consequences. 

This class contains 8 already obtained relations R^*" 30 ^ with polynomials (JSJ. 

Continuing the process of construction of decompositions and proper conse- 
quences we come finally to the prime relations R s *on<-2*3 defined over 4-simplices 
$i iii 2 i3 = {xi ,x llJ x l2J x l3 ,xg}, where i k <E {0, 1, ... ,7} and i < i x < i 2 < i 3 - 
The polynomials of these relations take the form 

P ' — XgO~4. (Xi , X-i-^ , Xi 2 , Xi 3 ) = XgXi Q Xi ± Xi 2 Xi 3 . (9) 

Substituting © in (©, |J5J, Q, and © (th is is a purely polynomial simplifi- 
cation) we have finally the following polynomial form of the system of relations 
valid for the Life rule: 



Xg {0-3 + a 2 + 1} + er 7 + 03 



(x 8 Xg + Xg) j a' 3 



3 ^ 
Xg {a\ 



0, 


(10) 


0, 


(11) 


0. 


(12) 


0, 


(13) 


0. 


(14) 
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Relations (|14|) have a simple interpretation: if the point xg is in the state 1, then 
at least one of any four points surrounding the center xg must be in the state 0. 

The above analysis of the relation R s Li f e takes < 1 sec on a 1.8GHz AMD 
Athlon notebook with 960Mb. 

To compute the Grobner basis we must add to polynomial © ten polynomials 



expressing the relation x p — x valid for all elements of any finite field F p ™ . 

Computation of the Grobner basis over F 2 with the help of Maple 9 gives 
the following. Computation for the pure lexicographic order with the variable 
ordering xg >- xs >- ••• >- xg remains initial polynomial @ unchanged, i.e., 
does not give any additional information. The pure lexicographic order with the 
variable ordering xo >- x\ >~ ■ ■ ■ >- xg gives relations ifTHjl — fH|) (modulo several 
polynomial reductions violating the symmetry of polynomials) . The computation 
takes 1 h 22 min. Computation for the degree-reverse-lexicographic order also 
gives relations (|10|l — (|14|) (with the above reservation). The times are 51 min 
for the variable ordering xg y x\ >- ■ • • >- xg, and 33 min for the ordering 

Xg y X 8 >~ ■ ■ ■ >~ Xq. 

4 Elementary Cellular Automata 

Simplest binary, nearest-neighbor, one-dimensional cellular automata were called 
elementary cellular automata by S. Wolfram, who has extensively studied their 
properties [7] . A large collection of results concerning these automata is presented 
in the Wolfram's online atlas [2]. In the exposition below we use Wolfram's 
notations and terminology. The elementary cellular automata are simpler than 
the Life, and we may pay more attention to the topological aspects of our 
approach. 

Local rules of the elementary cellular automata are defined on the 4-set 5 = 



function of the form s = f(p, g, r). There are totally 2 2 = 256 local rules, each 
of which can be indexed with an 8-bit binary number. 

Our computation with relations representing the local rules shows that the 
total number 256 of them is divided into 118 reducible and 138 irreducible re- 
lations. Only two of the irreducible relations appeared to be prime, namely, the 
rules 105 and 150 4 in Wolfram's numeration. 5 

We consider the elementary automata on a space-time lattice with inte- 
ger coordinates (x, t), i.e., x G Z = {...,—1,0,1,...} or x G 1 m (spatial m- 
periodicity), t G Z* = {0, 1, . . .} . We denote a state of the point on the lattice 

4 They are represented by the linear polynomial equations p + q + r + s + 1 = and 
p + q + r + s = for the rules 105 and 150, respectively. 

5 Wolfram prefers "big-endian" representation of binary numbers. 



x^ ~\~ x^, % — 0, . . . , 9, 



(15) 



{p, q, r, s} which can be pictured by the icon 
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by u(x,t) € S = {0, 1}. Generally the points are connected as is shown on the 
5x3 fragment of the lattice 




There are no horizontal ties due to the fundamental property of cellular automata 
— the states of points at a given temporal layer are independent. 

Applying our approach we see that some automata with reducible local re- 
lations can be decomposed into automata on disjoint unions of subcomplexes: 

1. Two automata and 255 are defined on disjoint union of vertices. 

2. Six automata 15, 51, 85, 170, 204 and 240 are, in fact, disjoint collections 
of zero-dimensional automata. What we call zero- dimensional automaton is 
spatially zero-dimensional analog of the Wolfram's elementary automaton, 
i.e., a single cell evolving with time. There are, obviously, four such automata 
with local relations represented by the bit tables 

1100, 

0110, (16) 

1001, 

0011. 



We call the automaton with bit table Hlfi[) oscillating point since its time 
evolution consists in periodic changing by 1 and vice versa. It is easy to 
"integrate" these automata. Their general solutions are respectively 

u(t) = 0, 

u(t) = u(0) + 1 mod 2, oscillating point, (17) 
u(t) = u(0), 
u{t) = 1. 

As an example consider the rule 15. The local relation is defined on the set 

® &) 

and its bit table is 0101010110101010. This relation is reduced to 

@> 

the relation on the face \j) and its bit table 0110 coincides with bit table 
(|16fl of the oscillating point. We see that the automaton 15 decomposes into 
the union of identical zero-dimensional automata on the disconnected lattice 




Using (|17fl we can write the general solution for the automaton 15 



u(x, t) — u(x — i, 0) + 1 mod 2. 

3. Ten automata 5, 10, 80, 90, 95, 160, 165, 175, 245, 250 are decomposed into 
two identical automata. 

As an example let us consider the rule 90. This automaton is distinguished 
as producing the fractal (of the topological dimension 1 and Hausdorff di- 
mension In3/ln2 w 1.58) known as the Sierpinski sieve, Sierpinski gasket, 

(P) (fft (r) 



or Sierpinski triangle. Its local relation on the set (j( is represented by 
the bit table 1010010101011010. The relation is reduced to the relation on 

® © 

the face (sf with the bit table 

10010110. (18) 

From the structure of the domain of the reduced relation it is clear that the 
lattice decomposes into two identical independent lattices as is shown 






To find a general solution of the automaton 90 it is convenient to transform 
bit table (|18fl to an algebraic relation. It is the linear relation s +p + r = 
and the general solution of the automaton takes the form 

u(x, t) = ^2 I J u(x - t + 2k, 0) mod i. 

k=0 ^ ' 

In the above examples we have considered the automata with reducible re- 
lations. If a local relation is irreducible but has proper consequences we also, in 
some cases, can obtain a useful information. 

For example, there are 64 automata 6 — both reducible and irreducible - 
having proper consequencies with the bit table 

1101 (19) 

on one or two or three of the following faces 



(20) 



The full list of these automata in the Wolfram's numeration is 2, 4, 8, 10, 16, 32, 
34, 40, 42, 48, 64, 72, 76, 80, 96, 112, 128, 130, 132, 136, 138, 140, 144, 160, 162, 
168, 171, 174-176, 186, 187, 190-192, 196, 200, 205, 206, 208, 220, 222-224, 234-239, 
241-254. 
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The algebraic forms of relation 119|) on faces (|20|) are 

ps + s = 0, qs + s = 0, rs + s = 0, 

respectively. 

Relation (|19|1 is non- functional. Nevertheless, it imposes a severe restriction 
on the behavior of the automata with such proper consequences. The peculiarities 
in the behavior are clear visible in the atlas [5], where many results of compu- 
tations with different initial conditions are pictured. A typical pattern from this 
atlas is reproduced in Fig.^ where several evolutions of the automaton 168 are 
presented. The local relation of the automaton 168 is pqr + qr + pr + s = 0. It 
has the proper consequence rs + s — 0. The black and white square cells in Fig. 
^ correspond to l's and 0's, respectively. Note also that the authors of Fig. ^ 
have used a spatially periodic condition. Their spacial variable is x G Z30. 











ipr ■ ■ ■ ■ f 





Fig. 1. Rule 168. Several random initial conditions 



Relation (|19fl means that if, say r, as for rule 168, is in the state 1 then s may 
be in both states or 1, but if the state of r is 0, then the state of s must be 0. 
Thus the corresponding diagonal or vertical may contain either only l's, or finite 
number of initial l's and then only 0's. The presence of a proper consequence of 
the form (|19H simplifies essentially computation with such automata: after the 
first appearance of 0, one can set 0's on all points along the corresponding line. 

In conclusion, let us present the results of analysis of the automata 30 and 
110. These automata are of special interest. The automaton 30 demonstrates 
chaotic behavior and even used as the random number generator in Mathe- 
matica. The automaton 110 is, like a Turing machine, universal, i.e., it is ca- 
pable of simulating any computational process, in particular, any other cellular 
automaton. The relations of both automata are irreducible but not prime. 

The relation of automaton 30 is 

1001010101101010 
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or in the algebraic form 

qr + s + r + q+ p = 0. 

It has two proper consequences: 

(p) @ (p) 

face 

bit table 11011110 11011110 

polynomial qs + pq + q rs + pr + r. 

The principal factor is 

1011111101111111 or qrs + pqr + rs + qs + pr + pq + s +p = 0. 

The Grobner basis of automaton 30 in the total degree and reverse lexicographic 
order is (omitting the trivial polynomials p 2 + p, q 2 + q, r 2 + r, s 2 + s) 

{qr + s + r + q + p, qs + pq + q, rs + pr + r} . 

We see that for the rule 30 the Grobner basis polynomials coincide with ours. 
The relation of automaton 110 is 

1100000100111110 (21) 

or in the polynomial form 

pqr + qr + s + r + q = 0. 

The relation has three proper consequences: 

(p) (g) (p) 

face 

bit table 11011111 11011111 10010111 

polynomial pqs + qs + pq + q prs + rs + pr + r qrs + s + r + q. 
The principal factor is 

1111111111111110 or pqrs = 0. 
The Grobner basis of automaton 110 contains different set of polynomials: 
{prs + rs + pr + r, qs + rs + r + q, qr + rs + s + q, pr + pq + ps} . 
The system of relations defined by the Grobner basis is: 

R {p,r,s} = 110lml = { prs + rs+pr + r = 0} , 
R {q,r,s} = 100lml = { qs + rs + r + q = Q y j 
R {q,r,s} = 10110 iii = { qr + rs + s + q = o} , 
R { P ,q,r,8} = 1110101110111110= {pr + pq + ps = 0}. 
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5 Conclusions 



Let us summarize the main novelties of the paper. 

— We have introduced a notion of a system of discrete relations on an abstract 
simplicial complex. Such a system can be interpreted as 

• a natural generalization of the notion of cellular automaton; 

• a set-theoretic analog of a system of polynomial equations. 

— After introducing appropriate definitions, we have developed and imple- 
mented algorithms for 

• compatibility analysis of a system of discrete relations; 

• constructing canonical decompositions of discrete relations. 

— We have proposed a regular way to impose topology on an arbitrary discrete 
relation via its canonical decomposition: identifying dimensions of the re- 
lation with points and irreducible components of the relation with maximal 
simplices, we define the structure of an abstract simplicial complex on the 
relation under consideration. 

— Applying the above technique to some cellular automata — a special case 
of systems of discrete relations — we have obtained some new results. Most 
interesting of them, in our opinion, is demonstration of how the presence 
of non-trivial proper consequences may determine the global behavior of an 
automaton. 

Acknowledgments 

This work was supported in part by the grants 04-01-00784 from the Russian 
Foundation for Basic Research and 2339.2003.2 from the Russian Ministry of 
Industry, Science and Technologies. 

References 

1. Toffoli, T.: Occam, Turing, von Neumann, Jaynes: How much can you get for 
how little? (A conceptual introduction to cellular automata). In: Automi Cellulari 
per la Ricerca e I'lndustria, Rende (CS), Italy, September 29-30, 1994- 

2. Hilton, P. J., Wylie, S.: Homology Theory: An Introduction to Algebraic Topology. 
Cambridge University Press, New York (1960) 

3. Lidl, R., Niederreiter, H.: Finite Fields. Encyclopedia of Math, and its Appl. 20, 
Cambridge University Press, London (1997) 

4. 't Hooft, G.: Dimensional Reduction in Quantum Gravity, Essay dedicated to Ab- 
dus Salam, Utrecht preprint THU-93/26, pp. 13; id., arXiv: gr-qc/9310026 

5. Bousso, R.: The Holographic Principle. Rev.Mod.Phys. 74 (2002) 825-874. arXiv: 
|hep-th/0203101| 

6. Gardner, M.: On Cellular Automata Self-reproduction, the Garden of Eden and 
the Game of 'Life'. Sci. Am., 224 (1971) 112-117 

7. Wolfram, S.: A New Kind of Scien ce. Wolfram Media, Inc (2002) 

8. |http://atlaa. wolfram.com/01/0l7] 



13 



